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^ Abstract 

I We investigate the production of a quark-antiquark pair in diffractive 

fl\^ photon-proton scattering, approximating soft pomeron exchange by the ex- 

r^ change of two nonperturbative gluons. In deep inelastic scattering at HERA, 

events with two jets and the scattered proton in the final state are predicted 
to be observable, with an important contribution from charm production. 
Eor photoproduction of light quark jets with high transverse momentum we 
find that both exchanged gluons must have a large invariant mass, so that 
the cross section is very small, whereas for charm quarks it is quite appre- 
ciable. Erom our calculation we also extract the quark structure function of 
the pomeron for the scaling variable z no too close to or 1, finding a strong 
fiavour dependence and a behaviour somewhat harder than z[l — z) for light 
quarks. 



1 Introduction 

Despite the great progress in describing strong interactions that has been 
achieved since the advent of QCD, processes that involve large distances, or 
small momentum transfers, are far from being well understood at present. 
Since the strong coupling becomes large at small momentum scales, pertur- 
bation theory cannot be directly applied in this situation, at least not in 
a straightforward manner. On the other hand, the phenomenology of soft 
processes, such as diffractive scattering or dissociation of hadrons, is fairly 
well described by the exchange of Regge trajectories. Most of these can be 
associated with mesons or baryons, but not so the trajectory with vacuum 
quantum numbers, the pomeron. The investigation of its properties and of 
how it can be described in terms of QCD is a challenge both for experiment 
and theory. 

At HERA, one is effectively colliding protons with real or virtual photons. 
One can look for collisions where the proton stays intact and loses only a tiny 
fraction of its momentum (instead of the proton there may also be a low- 
lying excitation such as Nit in the final state, cf. [1]), and where the pomeron 
it has radiated interacts with the photon. For purely kinematic reasons, 
the centre of mass of the photon-pomeron subreaction is well separated in 
rapidity from the scattered proton or proton remnant. Since the exchanged 
object is colourless there is no need for colour rearrangement in the final 
state which would "fill" this "rapidity gap" with further hadrons. Such a 
reaction can probe electromagnetic interactions of the pomeron. Given the 
high energies of the photon-proton collision at HERA, the photon-pomeron 
system can have an invariant mass large enough to allow for a hard reaction, 
which is amenable to a perturbative calculation, so that one is at the interface 
of soft and hard QCD. Events with a large rapidity gap and a jet structure 
in the final state have indeed been observed at HERA [2] and will be studied 
in more detail in the next few years. 

In this paper we investigate diffractive events with only the scattered 
proton and a pair of jets in the final state, the jets coming from a quark- 
antiquark pair produced in the photon-pomeron interaction (fig. 1). We 
consider both light and heavy fiavours for the quarks, and both photopro- 
duction and deep inelastic scattering (DIS). The framework of our calculation 
is the QCD-inspired model of Landshoff and Nachtmann [3, 4, 5], where the 
pomeron is approximated by two gluons with a nonperturbative propagator. 
It successfully reproduces some of the known phenomenological properties 
of the pomeron and has been used to describe exclusive diffractive p pro- 
duction in deep inelastic scattering. The process p + ^ -^ p -\- qq with real 
photons has previously been treated in this model [6], but there was an error 
in the calculation. Eor virtual photons the same reaction has been exten- 
sively studied in purely perturbative QCD, i. e. using the perturbative gluon 
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Figure 1: Diffractive production of qq by pomeron exchange. The figure 
corresponds to two diagrams because of the two possible directions of charge 
fiow in the upper quark fine. 



propagator [7, 8, 9, 10, 11]. 

In some aspects, the pomeron behaves quite hke a particle. Ingelman 
and Schlein proposed to describe processes in which a pomeron undergoes a 
hard reaction in terms of a parton structure of the pomeron and introduced 
the corresponding structure functions [12]. In the parton model for ordinary 
hadrons the quark structure function is related with the total cross section for 
photoabsorption, which can be generalised to the quark structure function of 
the pomeron and the cross section for diffractive photoabsorption. We can 
calculate the latter at Born level if we integrate the cross section for diffractive 
gg-production over the transverse quark momentum, and thus estimate the 
quark structure of the pomeron in our model. 

The present paper is organised as follows. In section 2 we describe the 
kinematics of the process under study and the kinematical limit we are work- 
ing in. Section 3 briefly recapitulates the Landshoff-Nachtmann (LN) model 
of the pomeron and section 4 describes the calculation of the cross section for 
diffractive gg-production in this model. Our results are presented separately 
for photoproduction and DIS in the following two sections. In particular, we 
flnd for photoproduction of light quarks that the exchanged gluons must each 
have a large invariant mass in order to contribute to the cross section. This 
is not the case for heavy flavours, and we obtain a cross section for photopro- 
duction of charm which should be observable at HERA (incidentally it is of 
the same order as the one calculated for light quarks in [6]). In section 7 we 
calculate and discuss the quark structure function of the pomeron in the LN 



model. The last section contains our conclusions. In an appendix, we briefly 
present what the phenomenological model of the pomeron by Donnachie and 
Landshoff [13, 14] predicts for our reaction, and show that it can be used to 
calculate the pomeron structure function, but not the production of high-px 
jets. 

2 Kinematics 

We begin with introducing some convenient kinematic variables for the pro- 
cess 

p{p) + i{q) -^ p{p) + q{P) + q{P') ■ (i) 

Let m be the mass of the quark q and v s the invariant mass of the quark- 
antiquark pair. We consider both real and virtual photons, and use the 
conventional abbreviations 

Q' = -q\ W' = {p + q)\ u = p.q (2) 

as well as Mp for the proton mass. 

For the momentum transfer of the proton (i. e. the momentum carried 
away by the pomeron) we use a Sudakov parametrisation 

p-p' = A = ^p + Cq + AT (3) 

with Ax transverse to both p and q. The on-shell conditions p^ = p'^ = M^ 
for the proton lead to a quadratic equation in ^ and ^' . Since we are interested 
in diffractive scattering, the squared momentum transfer t = A^ is small, say 
|t| < 1 GeV , and the ratio (W"^ + Q'^)/{s + Q^) is large, i. e. we are working 
in the Regge limit. As a consequence one finds that < (^ ^ 1, and that 
only one of the two solutions for the on-shell constraint is applicable, namely 

A|-2eMj 

t^^i-eMi. (4) 

Provided that s ^ |A^|, one can neglect !^' compared with !^. One then has 

and can interpret !^ as the proton's fractional loss of longitudinal momentum 
(cf. (3)). Requiring a minimum invariant jet mass ys translates into a lower 
bound for !^. We shall also impose an upper bound of 0.01 on !^ to be rea- 
sonably sure that the exchange of a pomeron dominates that of other Regge 
trajectories [6, 15], even at t sa 0. In most of our numerical applications, the 



cross section is dominated by smaller ^ and the effect of this bound is rather 
weak. 

In our calculation we shall neglect the transverse momentum At- To 
leading order in ^~^ one then has from (3) 

p '^ p '^ !^~^ IS. . (6) 

The proton mass will not appear in our results, since it turns out that terms 
involving Mp are suppressed at least by (^, such as in (4). 

We further introduce the transverse momentum pT of the quark relative 
to the photon direction in the qq centre of mass. Note that for fixed px 
and s the angle between the quark and the photon can take two different 
values which sum up to 180°. Interchanging these two values is tantamount 
to interchanging the quark and antiquark momenta, and due to charge con- 
jugation invariance the squared transition amplitude is the same for both 
configurations. Summing over them gives a factor of 2 in the phase space 
element. After integration over the azimuthal angles associated with pT and 
Ax, one has 

d^'p'd^^ifV' ,_ ^ ^^ 1 dtdjdpl . 

2p', 2Po 2P^ ^^ ^ y^ 2 3^i_4(p2^ + ^2)/j 

(7) 
We require the transverse jet momentum to be large, say |px| > 3 GeV 
at parton level, so that the jets can be identified experimentally and the 
produced quarks may be treated perturbatively in our calculation, including 
the case when Q^ = 0. 

3 The Landshoff-Nachtmann model 

In the LN model [3, 4, 5] the (soft) pomeron is modelled within QCD by 
the colour singlet part of two-gluon exchange. In contrast to a purely per- 
turbative approach, nonperturbative effects are partly taken into account 
by allowing for the confinement of the gluons or, more precisely, by work- 
ing with a nonperturbative gluon propagator. Nonperturbative effects are 
essential at low squared gluon momentum /^, whereas at hard scales pertur- 
bative behaviour dominates. When the perturbative piece, —g^^/l'^ in Feyn- 
man gauge, is subtracted at large P one remains with a gluon propagator 
—g^,yD[P) which can be related to the gluon condensate of the nonperturba- 
tive QCD vacuum. It was found in [3] that in order to reproduce the additive 
quark rule for total hadronic cross sections -D(O) has to be finite. Applied to 
certain processes, knowledge of the complete form of D[P) is not necessary. 



and only some of its moments are needed. In particular [4, 5, 6]^ 



'dl^[a<f^Di-F)] 
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dPP[a^^^Di-P)]' = ^^ , (8) 

Stt 

where /3o ~ 2.0 GeV~ is the coupling of the pomeron to u and (/-quarks [16], 
and ay> the strong coupling constant between on-shell quarks and a low-/^ 
gluon, which we will take to be ay> sa 1 (cf. [5]). The parameter /io ~ 1-1 GeV 
has been estimated from exclusive p production in deep inelastic scattering 
[4], and sets the characteristic scale over which D{1^) decreases. To the 
results obtained with this propagator one has to add the contribution from 
the perturbative piece which was subtracted at large gluon momenta. In 
diffractive hadron-hadron scattering it gives only a small correction [4]. 

Exchange of two noninteracting gluons is only the crudest approximation 
to the pomeron. We will calculate the leading term in ^~^ at one-loop level, 
and in order to make contact with experiment we shall modify by hand the 
exponent of !^ using the pomeron trajectory 

o/j.(t) = 1 -Ke-Ka't (9) 

with e ^ 0.08 [17] and a' ^ 0.25 GeV-^ [16]. 

The scattering of physical hadrons is calculated from the scattering of 
their constituent quarks. LN have found that, if the length /Xq is small 
compared with the hadron radius, the diagrams where the two gluons of 
the pomeron couple to different quarks in a hadron are suppressed, which 
naturally leads to the additive quark rule. Let us see what the remaining 
diagrams give for the radiation of a pomeron off a proton. 

We assume that the momentum of a constituent quark in the proton can 
be approximated hy p ^ xp with x ^ ^^ neglecting the transverse momentum 
of the constituents within the fast proton as we have previously neglected 
the transverse momentum Ax of the scattered proton. The product of Dirac 
matrices corresponding to the upper quark line in fig. 2 is 



t/3a 



^ + M,)7^(^-^+/+Mj7„(^+M,) , (10) 



where Mq is an effective mass for the constituent quark. Its precise value 
turns out to be irrelevant for our calculation, since corresponding terms are 
nonleading in ^~^ . Rewriting p k, xp ^ x^~^ ^ we see that to leading order 
inr' 

t^„ ^ (^+M,)7^(^+M,)7„(^+M,) 

^ xriA^(^+M,)7„(^+Mj + (af^/3) , (11) 



^The definitions in [4] difli"er from the ones used here by a colour factor of 2/9 (cf. next 
section), because the gluons are treated as abelian there. 
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Figure 2: A typical diagram contributing to the interaction between a proton 
and the pomeron, modelled by two nonperturbative gluons. The blobs stand 
for to the proton wave function and the dashed lines for gluons with Lorentz 
indices a and (3. 

where in the last step we have used the anticommutation relations for Dirac 
matrices. The matrix structure of this expression involves just the vector 
current of the proton, i.e. to leading order the spin structure of the pomeron- 
proton coupling is the same as for the coupling of the proton to a photon.^ 
Summing over the constituent quarks in the scattering amplitude one has for 
the proton vector current [18] 



{P\ E g(0)7«g(0)|p)^3Fi(t)u(/)7«^(p) 

q=u,d 



:i2i 



where the t-dependence is controlled by the isoscalar form factor of the nu- 
cleon, which we approximate by the Dirac form factor of the proton"^ 



Flit) 



AMI 



2.8t 



4M^ 



1 



t 



0.7 GeV^ 



(13) 



We will find in the next section that the factor x^ ^ in (11) is cancelled by 
the integration element of the loop momentum /, so that the coupling of the 



^Note, however, that due to its signature factor the pomeron has positive parity under 
charge conjugation, unlike the photon. 

^Since Fi(t) depends strongly on the squared momentum transfer t « A|, we do not 
neglect the transverse momentum Ay of the pomeron at this point. The above result (11) 
remains valid in this case if one replaces A with A — Ay 
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proton to the two gluons in the pomeron can be expressed in terms of A and 
the current in (12). 



4 Calculation of the cross section 

We shall now give some details on the calculation of the process p+7 — > p+qq 
or p + 7* -^ p -\- qq m the LN model. 

The leading term in ^~^ comes from the antihermitian part of the tran- 
sition matrix, which is equal to the imaginary part of the amplitude if an 
appropriate phase convention for spinors is used/ and can be calculated us- 
ing the unitarity relation for the scattering matrix. Fig. 3 shows the four 
corresponding Feynman diagrams. The lower fermion line stands for a con- 
stituent quark of the proton, which we shall treat as on-shell before and after 
the scattering. 

To project out the colour singlet part of the two-gluon amplitude, we note 
that the colour structure of all diagrams is given by 

rt^ (g) rt^ = -I (g) I - -t^ (g) t' , (14) 

t" denoting the Gell-Mann matrices divided by 2. Hence the appropriate 
factor for the colour-singlet amplitude is 2/9. 
For diagrams (c) and (0?) we parametrise 

l = lT + r]A + r]'q , (15) 

where 

is lightlike up to 0{^^) and It is transverse to A and q. Putting the four 
quarks coupled to the gluon with momentum / on shell gives a set of quadratic 
equations in r] and rj' . We introduce t = [q — Py and u = [q — P'Y as 
Mandelstam variables of the subreaction _P -|- 7 — > gg, where IP denotes the 
pomeron. The integration element then is, to leading order in (^~^, 

crt S[iP-iy-- ,n') . („ - A + ,y- - M^) = 5;^4^ • (17) 
One solution of the constraints gives, again to leading order. 



u — rw^ 



■*This result has been obtained in an approach where the exchanged gluons are treated 
perturbatively [8]. For a nonperturbative gluon propagator it can be reproduced by a 
calculation along the lines of [3], where only analyticity properties of D{P) are needed. 
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Figure 3: Feynman diagrams contributing to the imaginary part of the amph- 
tude for p+7 — > p+qq. The pomeron is approximated by two nonperturbative 
gluons (dashed hues). 



whereas the other leads to 

/2 = x^-^iu-m^) , (19) 

which is very large due to the factor ^~^ and therefore has to be discarded. 
For the other two diagrams we make a similar ansatz 

i = -It + v^ + v'q (20) 

and, choosing the sign of It in this particular way, obtain analogous results 
with ii replaced by t and P by P' in equations (17) to (19). For the denomi- 
nator of the off-shell quark propagator we then find da = t — m^ , dc = ii — m^ 
and 



dh dd 1 



u — m'^ t — rrr Py -|- m^ s -\- 



(It - PT? + ? (P^ + m') + m' 
s 



(21) 
for diagrams (a), (c), (6) and (o?), respectively. The kinematics of the first 
two of these diagrams are clearly different from those of the other two, and 
the contributions from the two separate sets will be discussed is some detail 
below. 

Together with (11) and (12) for the coupling of the gluons to the proton 
it is now straightforward to calculate the differential cross section for our 
process.^ We shall give some intermediate results which will be useful in the 
discussion. 

Let a be the electromagnetic fine structure constant and Cq the charge of 
the produced quarks in units of the elementary charge. In each diagram of 
fig. 3 there is one quark-gluon vertex which involves an off-shell quark. We 
take the running strong coupling as{pj) instead of a^' for this vertex in all 
four diagrams, with the usual uncertainty about which argument should be 
chosen for the running coupling to keep higher loop corrections small. In our 
numerical studies we use a fixed Q^gij)^) = 0.2, given that the variation of 
the running coupling is not very large in the relevant range of \)\. We use a 
fiux factor 2{W'^ -\- Q^) = 2!^~^{s -\- Q^) for the photon, take the average over 
the spin of the incoming proton, and sum over the spins of the final state 
fermions as well as the colour of the produced qq pair. Then 



d(T 



T,L 



dt d^ dp^ 
16a,s{pT 



a 



<°'-'^'"°""'^""'lW Vl-4(p| + ^-)/^ ^--''"' 



®We have used the algebraic programs FORM by J Vermaseren and MATHEMATICA 
by Wolfram Research, Inc. 
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where 



Xt 



+ 



p| + m^~ 




1 



Pj, + m 



Pj + m^ 



- [«i°)i^(-l^)]^ 



(/2/, 



Pt + (It - Pt 



TT 



u — m 



u — m 



dh 



dh 



(23) 



and 



It 



^2 ^2 



p^ + m^ If d^l 



s + 



2g^ 




- [«i°^i^(-i^)]^ 



ii — m^ 5 + 0^ 



0?h 



2g2 



(24) 



for photons with transverse or longitudinal polarisation in the proton-photon 
centre of mass, respectively. The cross section is the same for the two possible 
transverse polarisations because we have integrated over the azimuthal angle 
of the jets. 

Note that it is possible to write the longitudinal cross section as the square 
of an integral over one loop momentum.^ Similarly, the transverse cross 
section is the sum of two squared integrals, the integrand in the first term 
being a vector quantity. By calculating separately the squared amplitudes 
for the sums of diagrams (a), (c) and (6), (o?), as well as their interference 
term, it is possible to separate the contributions of these two sets to the 
cross section. In (23) the first term in each of the large brackets comes from 
the sum of (a) and (c), the second (with o?6 in the denominator) from the 
sum of diagrams (6) and (o?). Likewise, in (24), the first expression enclosed 
in curly braces comes from (a) and (c), the second one from (6) and (o?). 
For the latter, the denominator o?6 has cancelled in the term proportional to 
l/Q^- One can see that the two sets of diagrams are not separately invariant 
under the electromagnetic gauge group, as the amplitudes for the sums of 
(a) and (c) and of (6) and (o?) diverge for Q^ — > 0. Similarly, it turns out 
that these sums are nonzero if the photon polarisation is set proportional to 
its momentum q. 

One can easily perform the integrations over the azimuthal angle of It 
and then has 



®It cannot readily be identified as proportional to the amplitude because the cross 
section is summed over polarisations and azimuthal angles in the final state. 
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Pj + rw^ \Jo J 

U = 4^P^:+^(rjiM«(°)i)(-i^)]V2(^,-))' , (25) 

s s \Jo J 

where we have introduced the functions 

u + ti; — 2 



w 
h{v,w) = 1 + — 



^{v + w-2y -4(1 -w) 



w 



h{v,w) = 1- I (26) 

^J{v + w -2f -4(1 -w) 



of the arguments 

1| pl + m^ s + Q^ 

y = —- ^ w = — . (27) 

Pt Pt ■s 

Fig. 4 shows fi[v^w) and f2{v^w) for some fixed values of w. Independent 
of w^ both /i and /2 vanish at u = and tend to 1 as u — > oo. For w = 1^ 
/i becomes the step function 6[v — 1). The singularity of /2 at u = 1 in the 
same limit is spurious, since w only becomes 1 if both Q^ and m^ go to zero, 
and the integrals with /2 are multiplied with one of these factors. 

This representation of the cross section is similar to the one obtained by 
Nikolaev and Zakharov [9], who work with a perturbative gluon propagator. 
Our equations (22) and (25) to (27), taken at t = 0, are in agreement with 
eq. (22) to (27) of ref. [9]^ if 

1. we replace ap[t) in our eq. (22) with 1, i. e. if we do not introduce the 
pomeron trajectory in order to account for the interaction between the 
gluons, 

2. one takes into account that the strong coupling is taken at different 
momentum scales in the present paper and in [9], 

3. our [_D(— 1|.)]^ is replaced with V{\t)/(\-t + /"g)^5 corresponding to the 
different assumptions in [9] about the gluon propagator and the cou- 
pling between the gluons in the pomeron to a hadron. 



"^We assume that in the second line of their eq. (24) it should read k^ instead of k. 
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(b) 



Figure 4: The functions /i and /2 of eq. (26) for w = 1.1 (full), w = 2 
(dashed) and w = 10 (dot-dashed). /2 has a zero at u = 2(2 — w). 
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4. one makes a transformation from our integration variable ^ to their a 
with 

-5 = 4±^ (28) 

a[L — a) 

(eq. (19) of [9]). Using the relation (5) between ^ and s we find 

da = 2ds ^-^^ 1-4 ^-^ 

s^ \ s ) 

= 2t/^r ^^ 1-4^*-^^^ . 29 

s s \ -^ J 

Note that for fixed \)\ and s there are two values of a, which gives 
a factor of 2 in the integration element as explained at the end of 
section 2.® 

Assuming that due to the squared gluon propagator [_D(— 1^)]^ the main 
contribution to the loop integrals comes from values of \\ up to the order of 
/Xq sa 1.2 GeV , which is small compared to the minimum \)\ required for the 
jets, we can expand /i and /2 to first order in ly/py and obtain, using the 
second moment in (8) 



n 
?4,,4 P.2 



ll 



81/3oVo s^ 1 




(s + g2)2 (p2^ + ^2^2 

p^ + m 



m 



s 

2 



4 Pj + m'^ 

16^2 (J + 




^L - ^^7S^7-^^7^Z^^ri(l-2^^Z2^^^) • (30) 



To assess how good this approximation is, we have evaluated the loop inte- 
grals in (25) for a specific model of the gluon propagator, namely [4] 

where n > 4 and the proportionality constant can easily be obtained from 
(8). For n — > oo this becomes D[—\^) oc exp(— l|.//iQ). The approximation in 
(30) improves slightly with increasing n. Fig. 5 shows that it is rather good, 
although at small Q'^/s it underestimates the cross section for transverse 
photons.^ 



®We differ from eq. (6) of [11] by this factor. 

^If at Q^ = 5 GeV , s = 400 GeV and m = 0, we integrate the transverse differential 
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2 \n^\r'i] 



Q" [GeV 




100 



g2 [GeV2 



(&) 



Figure 5: (a) The squared loop integral of eq. (25) which involves /i (full), 
compared with its approximation in (30) (dashed) for |px| = 3 GeV, s = 4p|., 
m = and n = 4. To obtain the correct Q^-dependence of the cross section 
it has been multiplied with s/[s + Q^). Up to this global factor, changing s 
and m corresponds to rescaling Q^ since /i depends on these parameters only 
through w in (27). (6) Same as (a) but for /2 instead of /i and with a global 
factor of Q"^ I {s + Q^) instead of s/(s + Q^), so as to obtain the Q^-behaviour 
of the longitudinal cross section. 
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In the calculation of ref. [6] the propagator denominators dj, and 0?^, which 
depend on the loop momentum, have been approximated by their values at 
It = 0. This leads to a wrong result, because at 1^ = the integrands vanish, 
and to obtain the terms of order \'^ the propagators have to be expanded to 
this order, too. 



5 Results for photoproduction 

We will now discuss our results for the case of photoproduction, concentrating 
on the production of light quarks. Setting m = m eq. (30) one finds that the 
cross section for transverse photons is proportional to Q"*, and thus vanishes 
for photoproduction even faster than the longitudinal one which behaves like 
Q^. This is not an artifact of the linear approximation in l^/px which led 
to (30). At m = and Q^ = the function /i in (25) is the step function 
6{\^ — pf^) and the loop integral is strongly suppressed by a factor of order 
[-D(— p^)]^. For squared momenta of order p^ the perturbative part of the 
gluon propagator is of course dominant. To calculate its contribution one 
can substitute l/{—\j) for D{—\j) and as{pj) for ay\ taking all four quark- 
gluon vertices in fig. 3 as perturbative and using the perturbative coupling 
at the fixed scale p'^ for simplicity. Then the cross section for m = and 
g2 = becomes (cf. (22) and (25)) 

dcr 



dt d^ dp J' 



"^ V 1 ~ 4p2-/s \ ■s / Pt 



With |px| > 3 GeV and ^ < 0.01 this gives a very small integrated ^p cross 
section of 9.5 pb for u and o?-quark jets at a photon-proton centre-of-mass 
energy of W = 200 GeV. The effect of imposing an upper limit in the inte- 
gration over ^ is of the order of 10%. 

We must however bear in mind that because of the highly virtual gluons 
we are no longer dealing with the soft pomeron known from phenomenology. 
The introduction of the soft pomeron trajectory ap[t) to take into account 
higher order corrections to the exchange of two noninteracting gluons in fig. 3 
is presumably no longer justified. We will not attempt to discuss here the 
difficult problem of the interplay between the soft pomeron and its pertur- 
bative counterpart, but only mention that if one substitutes for example the 
term ^~'^^ in (32) which comes from the soft pomeron trajectory with ^~^ 
corresponding to the Lipatov value for the pomeron intercept, the integrated 



cross section over p|, from 9 GeV to its upper kinematical limit, the approximation gives 
70% of the exact value for n = A. 

16 



cross section mentioned above becomes as large as 1.8nb.^° Another question 
is whether the couphng of the Lipatov pomeron to the proton is dominated 
by soft or hard scales. In the approximation of two noninteracting gluons, 
the quark-gluon coupling has of course to be perturbative at the bottom of 
the diagrams in fig. 3 if it is perturbative at the top. 

Let us analyse in some detail what leads to the suppression of loop mo- 
menta \j < pj for photoproduction and massless quarks. Eq. (21) for the 
denominators of the off-shell quark propagators in diagrams (6) and (o?) now 

reads 

ch ^ck ^ (It - PtY .gg. 

u i Pt ' 

and from eq. (23) we take that the loop integral to be calculated is 



IT 



[d{-\'t)Y 



p2 

Pt + (It - Pt) ^ 



(It - Pt^ 



(34) 



or the analogous expression with 1/1^ instead of [_D(— 1^)]^. The first term 
in large brackets comes from diagrams (a) and (c), the second from (6) and 
(d). The integral for the contribution from (6) and (d) is well known:^^ for 
It ^ Pt i^ averages to zero when the azimuthal angle is integrated out, so 
that only (a) and (c) contribute to the cross section, whereas for \'^ < p'^ it 
gives a finite contribution which just cancels that from (a) and (c). 

Since the squared gluon propagator decreases rapidly with 1^ the loop 
integral is dominated by the region where \j is not too large compared with 
p'^. In a small part of this region, namely when It ~ Pt, it follows from 
(33) that in diagrams (6) and (o?) the real photon decays into two quarks 
which are on or near shell and as a consequence have to be approximately 
coUinear with the photon direction. In this case one can expect important 
strong corrections to the coupling between the photon and the quarks, which 
in our calculation was taken as the electromagnetic Born-level vertex. 

Let us finally investigate the effect of a nonzero quark mass. For light 
quarks, the contribution to the cross section from the term with /i in (25) is 
completely negligible for the same reason as above, since the parameter w in 
(27) is very close to the value 1, which it takes in the case of zero mass. There 
is, however, a small contribution from the term proportional to m^ which in- 
volves /2. Taking m = 150 MeV for strange quarks the nonperturbative con- 
tribution to the integrated cross section is very small (sa 2pb). As indicated 
at the beginning of this section, the transverse cross section calculated with 
the nonperturbative propagator D[P) vanishes faster than the longitudinal 



^"^Here we have used the slope a' of the soft pomeron, but the numerical effect of 
changing its value is not too great. 

^^In two-dimensional electrostatics it occurs e. g. when the electric field of a rotationally 
symmetric charge distribution is calculated from Gauss' law. 
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one for Q^ — > 0. For production of u, d and s quarks at Q^ = 1 GeV the 
former is 3.4 pb whereas the latter is 45 pb, with W = 200 GeV, |pt| > 3 GeV 
and ^ < 0.01. 

To see how sensitive our results are to the choice of taking current masses 
for the produced quarks we have also done the calculation with constituent 
masses, 300 MeV for u and d and 450 MeV for s. Then, due to the term pro- 
portional to m^ in (25), the nonperturbative contribution to photoproduction 
of the three light quark species is 51 pb with W, |px| and ^ as above. 

In the case of heavy quarks the situation changes and becomes similar to 
that in DIS. Photoproduction of charm quarks will therefore be treated in 
the next section. 

6 Results for deep inelastic scattering 

When Q^ is sufficiently large, small \'^ are no longer suppressed in /i, because 
w in (27) is larger than 1. The same is true even at low Q^ for a large quark 
mass, such as for charm quarks with m = 1.5 GeV. From 



i-m' = -'^^Y~ (l ± ^1 - 4(Pt + m^)l-^ (35) 

and the corresponding equation for u with the opposite sign in front of the 
square root, one can see that the denominators (21) of the propagators in 
diagrams (6) and (o?) are limited by 



\d,l\dd\>Q' ^''\'^\ m' . (36) 



Therefore the corresponding quarks are well off-shell and the use of the Born 
level coupling with the photon, as well as that of the perturbative quark 
propagator, is justified. 

Let us us now discuss some numerical features of the cross section cal- 
culated in the LN model. Table 1 gives integrated cross sections obtained 
from our approximate formulae (30). For light quarks, the longitudinal cross 
section (Tl is larger than the transverse one at very low Q^ (which is still 
an effect of the low-l|. suppression in /i), but the transverse contribution 
overtakes between Q^ = 10 GeV and Q^ = 15 GeV . The transverse cross 
section a^ reaches its largest values around Q^ = 40 GeV , and becomes 
quite small beyond 100 GeV^. For longitudinal photons, the cross section 
decreases for Q^ > 10 GeV and is rather small when a^ takes its maximum. 
We should remark that the cross sections given in table 1 are lowest-order 
results which will be modified by perturbative QCD corrections to the upper 
part of the diagrams in fig. 3, in particular at large values of Q^. 
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Table 1: Integrated photon-proton cross sections obtained from (30) with a 
transverse jet momentum |px| > 3 GeV and a fixed photon-proton centre- 
of-mass energy o{ W = 200 GeV. The effect of the cut !^ < 0.01 is rather 
smaU. 



Whether one takes current or constituent masses for the fight quarks does 
not make mucfi difference for tfie resufis. Tfie mass dependence is largest 
at low Q^, wfiere taking constituent masses one obtains a somewfiat larger 
transverse cross section. For Q^ > 15 GeV tfie cfiange in a^ wfien taking 
current or constituent masses is below 10%, and tfie same is true for tfie 
different values of (Tt for u, o? or s quarks (apart from tfie different electric 
cfiarges). In (Tl tfie effect is yet smaller at all values of Q^. 

For cfiarm quarks, tfie results are quite different. As we fiave mentioned in 
tfie previous section, tfie transverse cross section does not vanisfi at Q^ = 0. 
In fact, it takes its largest value tfiere and decreases witfi Q^ at fixed W. 
Tfie longitudinal cross section is negligible for all values of Q^. 

As for tfie dependence on tfie centre-of-mass energy, tfie cross section 

rises slowly witfi W. As an example, witfi Q^ = 50GeV^, |pt| > 3 GeV and 

^ < 0.01 we obtain a transverse cross section for u, d and s production (taking 

current masses) of 210 pb, 280 pb and 300 pb for W = 140 GeV, 200 GeV and 

250 GeV, respectively. Tfiis befiaviour can be understood if, using (5), one 

makes a transformation 

d^ ds 



(37) 
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and rewrites (22) and (25) in the form 

da OC ^2(l-«r(t)) [^^(^)]2 J^ JJ Jp2^ ^ (3g) 

where the constant of proportionahty depends on the variables s, p'^ and Q^ 
only. Expressing ^ in terms of W^ and s, we see that the only dependence on 
W comes from the term |t2(i-«2p(*)) with its small exponent. The biggest effect 
is actually that the phase space is reduced at small W. We impose an upper 
limit on (^, and as a consequence there is a maximum allowed s increasing 
with W. If we integrate without any upper bound for (^, the above cross 
sections become 270 pb, 290 pb and 300 pb with an even weaker dependence 
on W. We also see that except for rather low values of W the effect of the 
cut on ^ is small, for W = 200 GeV the results with or without cut differ by 
less than 4%. 

Some features of partially integrated cross section are shown in figures 6 
and 7. The (^-spectra in fig. 6 show that for W = 200 GeV the cross section is 
quite small at ^ = O.OI. One might remark that the differential cross section 
for longitudinal photons has a zero for certain values of the parameters. 
This can be most easily seen from the term in large brackets in (30), which 
vanishes at Q^ = s for m = and at smaller Q^ if m ^ 0. Eq. (30) is only an 
approximation, and the true position of the zero is slightly different (fig. 5). 
That a zero does exist also follows from the fact that /2 changes its sign at 
low \'^ as w increases (cf. fig. 4), so the loop integral involving /2 changes 
sign as well and must have a zero for some value of w. Due to the rather 
high values of Q^ at which this occurs it will however be very difficult to see 
the zero in experiment, because the transverse cross section is comparatively 
large in this region. 

The p|.-spectra in fig. 7 are fairly well described by a power behaviour 
da I d\)\ OC |px|~ . Eor production of the three light quark fiavours with trans- 
verse photons, average values of the exponent are 8 sa 8.7 at Q^ = lOGeV^ 
and S sa 6.4 at Q^ = 40 GeV , and with longitudinal photons we find S ^ Q at 
Q^ = 10 GeV . The spectra for charm quarks are less steep, with exponents 
5 ^ 6.7 at g2 = lOGeV^ and 5 ^ 5.7 at Q^ = 40 GeV^ both for transverse 
photon polarisation. 

We finally estimate how the contribution from the perturbative gluon 
propagator at large \'^ will modify these results. There is some arbitrariness 
in choosing the value of |lx| above which it should be taken into account [4], 
but it should be at the GeV-scale, and thus of the same order as the mini- 
mum |px| required for the jets. Using a lower cutoff for the 1^-integration of 
k"^ = 4GeV^ and replacing a(°)_D(— 1|) with as{\j)/{—\j),^^ we have stud- 
ied the ratio between the perturbative and the nonperturbative contribu- 
tions to the differential cross section (22), (30). In the case of transverse 



^^One can also take the running coupling at p|, instead of 1|, for the vertex involving 
the off-shell quark. The numerical difference is very small. 
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Figure 6: (a) dcrx/d^ summed over u and d quarks [m = 0, |pt| > 3 GeV, 
W = 200 GeV) at Q^ = 20 GeV^ (full), 50 GeV^ (dashed), 100 GeV^ (dot- 
dashed), (b) The same for da^/d^ at Q^ = 5 GeV (full), 10 GeV (dashed), 
20GeV^ (dot-dashed). 



21 



do-/ dp J. [pb GeV 

100 




p| [GeV2 



da /dp J. [pb GeV 

100 




p| [GeV^ 



Figure 7: (a) The p|-spectrum da/dpl for ^ < 0.01, W = 200 GeV and 
Q^ = 10 GeV . The curves are for the production of u, d and s quarks (with 
current masses) and transverse 7* (full), the same quarks and longitudinal 
7* (dash-dotted), and for charm quarks and transverse 7* (dashed). (6) The 
same at Q^ = 40 GeV , without the curve for the (rather small) longitudinal 
cross section. 
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photons and light quarks it increases with p'^ and is largest for low Q'^/s^ 
where the nonperturbative cross section is very small. If one sets m = 
and s = 4p^, then for |px| = 3 GeV the perturbative tail gives a correc- 
tion of 10% for g2 = 5GeV2 and less than 3% when Q^ > 15GeVl With 
|Pt| = 6 GeV one finds an 8% correction at Q^ = 40 GeV , which grows to 
19% at Q^ = 15 GeV . There however, the longitudinal cross section domi- 
nates over the transverse one anyway, so the perturbative correction is not 
too large at least for values of Q^ where the transverse cross section is impor- 
tant. Also, the absolute size of the perturbative contribution is always small 
compared with the scale set by the maximum values of the nonperturbative 
transverse and longitudinal cross sections. As for the production of charm 
quarks, we have seen that the nonperturbative transverse cross section does 
not vanish for low Q^. Even at Q^ = the perturbative correction is below 
4% for |px| = 3 GeV, and less than 11% when |px| is twice as large. 

In the case of longitudinal photon polarisation the perturbative correc- 
tion to the cross section is smaller and can be neglected. Where the cross 
section for longitudinal photons calculated with ay>D{—\j) vanishes, the 
contribution from the perturbative tail is unobservably small compared with 
the transverse cross section. 

As remarked in the previous section, the simple substitution of the non- 
perturbative by the perturbative gluon propagator and the introduction of 
the running strong coupling for all vertices may not be sufficient to describe 
the physics at large gluon momenta. If one takes the Lipatov intercept of the 
pomeron trajectory when the gluons in our one-loop calculation are hard, the 
contribution from the perturbative tail dominates over the one from small 
\'^ when ^ is some parts in 10""^ in the case of transverse photons, for lon- 
gitudinal ones it is larger for some values of the parameters p|., s, Q^, and 
smaller for others. We emphasise again that this is not meant as a proper 
treatment of the perturbative pomeron, but rather as an indication of what 
the numerical situation is under a certain assumption. 

7 The pomeron structure function 

In the present section we shall use our calculation of diffractive jet production 
to estimate the quark structure function of the pomeron. 

More precisely, one can generalise the usual relation between the cross 
section for the absorption of transverse photons by the proton and the quark 
distribution functions qi^xs), 

-T{l*p) = ^^,Y.elq{xs) (39) 

with our choice of 2{W'^ -\- Q^) for the photon fiux factor, to the differ- 
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ential cross section dcTT I {dt d!^) for diffractive photoabsorption, introducing 
the diffractive distribution functions dq/[dtd^). Here xb = Q'^/{2i') is the 
Bjorken variable, and the sum over q is to be taken for quarks and anti- 
quarks. Under the assumption that diffractive scattering is dominated by 
pomeron exchange, and that the couphng of the pomeron to the proton and 
the diffractive final state factorises, one can write dq/[dt d^) as the product of 
the "fiux of pomerons in the proton" and a quark structure function Gq/p{z) 
of the pomeron. Since the pomeron is an isosinglet, the quark and antiquark 
distributions are equal. We use the definition of [f3], that is 

where 

z = xb/^ = QV{s + Q') (41) 

is the scaling variable for the photon-pomeron subreaction. Together with 
the diffractive equivalent of (39) this gives 



dar ex 



9(3nFi{t)re-'"^^'^j:4-GM^) ■ (42) 



dtdc s + g2 ' UL .V /J ^ , ' z 



In the phenomenological model of the pomeron which is described briefly in 
the appendix, one obtains a structure function 

G,/Az) = ^^ ^(1 - ^) ~ 0.2 zil - z) (43) 

for each massless quark or antiquark and z not too small, say for 2; > 0.1 

[13]. _ 

Since p + 7 — > p + gg (flg. 1) is the "Born approximation" of diffractive 
photoabsorption, we may use our present calculation to estimate the pomeron 
quark structure in the LN model. This requires integration of the differential 
cross section (22), (25) over p^ from zero to its kinematical upper limit. A 
hard scale in the photon-pomeron reaction is now provided by Q^, and we 
shall furthermore require that the squared invariant mass s be large enough, 
say s > 20 GeV . For flxed Q^ this puts an upper limit on the scaling variable 
z in (41) for which we assume our calculation to be sensible. However, it can 
be seen from (21) and (35) that if p^ is small and Ix is near px, the squared 
momentum of the off-shell quark in diagram (6) or (0?) becomes small for light 
quarks, so that the perturbative treatment of the quark-photon coupling 
and the quark propagator in the upper part of the concerned diagram is 
problematic in a part of the region of integration. 

We have calculated Gq/p{z) in 2;-steps of 0.1 with Q^ equal to 5, 10, 20, 50 
and 100GeV^ the requirement s > 20 GeV^ leading to z < 0.2, 0.3, 0.5, 0.7 
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and 0.8, respectively. We further take z > 0.1, corresponding to ti; > 10/9, 
so that the problems associated with w ^ \ which occur for photoproduction 
do not exist (cf. fig. 4).^"^ 

When calculating the p|.-integrated cross section, we can no longer ex- 
pand the loop integrand in ly/py, after which only a moment of the gluon 
propagator was needed, and therefore we take the particular form (31) of 
D[P) for different values of n. The dependence of our results on n is very 
weak, of the order of 5%, which gives some hope that the precise shape of 
D[P) is not too important. As for the running coupling as(p|.), which we 
introduced for the vertex involving the off-shell quark, we choose to freeze it 
when it becomes equal to one, admitting that again this is far from being a 
rigourous treatment. Using the QCD scale parameter A^ = 260 MeV [19], 
the freezing point is at 570 MeV and varies by about 100 MeV when the cited 
errors on A^ are taken into account. 

We find that our results are not very different in size from (43), with 
a shape that is somewhat shifted towards larger z. In the following, we 
shall give the ratio of Gq/p{z) obtained in the LN model and Gq/p{z) in 
(43). Fig. 8 shows this ratio for light quarks. The scaling in z is very good: 
where we have several Q^ for the same z, differences in the structure function 
are below 1%. Looking at (22) and (25) to (27), we observe that the loop 
integrals depend on s and Q^ only via w^ and thus directly on z^ and that 
the additional s-dependence via (p|. -|- m^) j s is not very important for small 
values of p^. In fact, low p^ dominate in the integrated transverse cross 
section, the contribution from the region p^ > /Xq being less than 10% for 
?7i = and less than 20% for m = 450 MeV. The dominance of low p|., 
which is already seen in jet production (cf. fig. 7), can be understood from 
the behaviour of /i, which is increasing in u = \\l\)\. For /2 the situation 
is less clear-cut, but the corresponding term in the transverse cross section 
is suppressed by m^ j {j)\ -\- m^) and accounts for at most 25% of the result 
for light quarks. With the integrated cross section being dominated by low 
p|., its size is infiuenced by the value we assume for the strong coupling a^' 
rather than by the scale from which onwards we take the running coupling 
o^s(Pt)- If foi" example one uses a^ = 1 for all values of p^, the result for the 
structure function grows only by a factor between 1.2 and 1.4 at ?7i = and 
by 1.4 to 1.6 for m = 450 MeV. 

Because we require s not to be too small in our calculation, we need to go 
to rather high Q^ to cover the whole z range (cf. (41)). The scaling we find 
presumably exists only in a lowest order approximation; it is to be expected 
that there will be scaling violation as Q^ increases, due to loop corrections 
to the upper parts of the Feynman diagrams (fig. 3) or, in other words, due 



^^Also, for very small z, the contribution from the triple pomeron vertex becomes im- 
portant [13]. 
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Figure 8: The ratio of the pomeron quark structure function Gq/p{z) ob- 
tained in the LN model and the result (43) from the phenomenological ap- 
proach of ref . [13]. The points correspond to Q^ = 100 GeV and 171 = (o), 
m = 150 MeV (+), m = 300 MeV (n), and m = 450 MeV (x), respectively. 
Points for other values of Q^ coincide with the ones shown here to a good 
precision. 

to the perturbative evolution of the structure function. 

The dependence on the quark mass is rather pronounced and leads to 
a flavour asymmetry. The structure function is smaller for s quarks than 
for u and d quarks, and there is a difference both in shape and in the nor- 
malisation according to whether one takes current or constituent masses. 
Increasing m leads to a suppression at low p'^ of the /i-term in (25), due 
to the factor Px/(px + rn'^)^ and at the same time increases the parameter 
w oc (p|. -|- ?7i^)/p|., SO that /i itself becomes smaller for most values of v 
(cf.fig. 4). 

From dcrL/[dt d^) one can obtain a longitudinal structure function anal- 
ogous to (42), which for light quarks turns out to be relatively unimportant. 
The longitudinal cross section is suppressed at low p'^ by a factor [p'^-\-m'^)/s^ 
and the longitudinal structure function is less than 6% of the transverse one 
for m = 0. Its importance grows with increasing m: at m = 450 MeV the 
ratio of longitudinal and transverse cross section reaches 20% for large z^ but 
is still under 6% if one takes out the points at Q^ = 50 GeV , z = 0.7 and at 



Q^ = 100 GeV^ 



0.8. 



For charm quarks, Gq/p{z) is about an order of magnitude smaller and 
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Figure 9: The pomeron structure function for charm quarks divided by (43). 
Displayed is the transverse structure function at Q^ = 100 GeV (□), which is 
almost the same for other values of Q^. The longitudinal structure function 
is shown at Q^ = 50 GeV^ (+) and at Q^ = 100 GeV^ (x). 

suppressed at large z (fig. 9), for the same reasons discussed above in connec- 
tion with light quark masses. The dominance of low p^ is less pronounced, 
values of |px| above the charm mass contributing between 12% and 35% to 
the integrated cross section. As a result, the perturbative region of a^ is now 
more important (taking a constant a^ = 1 gives an answer between 2 and 2.5 
times as large). The scaling in z is still quite good: the variation of Gq/p{z) 
with Q^ is at most 12%, and less than 7% when the values at Q^ = 5 GeV 
are not taken into account. The importance of the longitudinal structure 
function is no longer negligible at large z. It decreases with Q^ at fixed z^ 
i. e. it shows no scaling, and increases strongly with z. For comparison the 
longitudinal structure function at two different Q^, normalised to Gq/p{z)^ 
is included in fig. 9. 

Finally we remark that the contribution of the perturbative tail of the 
gluon propagator to the structure function can be neglected. The latter is 
dominated by low p'^ and, as we have remarked at the end of section 6, the 
importance of the perturbative corrections is most important at high p'^. 
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8 Conclusions 

In the present paper we have calculated the cross section for diffractive pro- 
duction of a quark-antiquark pair with high transverse momentum px in 
photon-proton scattering, modelling the pomeron by two noninteracting glu- 
ons with a nonperturbative propagator D[P). To a good approximation the 
result only depends on a moment of D[P) which has been estimated from 
exclusive diffractive p production. 

We find that in deep inelastic scattering at HERA such events should 
occur with an observable cross section. The rate for charm quark produc- 
tion is comparable to the one for the light quark fiavours, a prediction that 
appears to us worthwhile testing in experiment. For light quark production 
there is an important contribution from longitudinally polarised photons at 
moderate Q^, up to about 20 GeV . As expected for soft pomeron exchange, 
the dependence of the cross section on the photon-proton centre-of-mass en- 
ergy W is rather weak, the largest effect coming from the reduction of phase 
space at low W, which is due to an upper cutoff we impose on the fractional 
momentum loss ^ of the scattered proton. 

In the case of diffractive photoproduction of light quark jets, there is a 
cancellation in the loop integration unless the gluons are hard, at least if 
one takes current masses for the produced quarks. To give a nonvanishing 
contribution, the squared invariant mass of the gluons, which is equal to their 
squared transverse momentum relative to the pomeron direction, has to be 
at least of the order of the squared transverse momentum of the jets, and 
as a consequence the nonperturbative part D[P) of the gluon propagator 
gives an unobservably small cross section. The diagrams we calculate give 
only the lowest order approximation of pomeron exchange in the process 
under study. We do not know how corrections from higher loops will modify 
this result, in particular the replacement of two noninteracting gluons by 
a ladder structure, which is believed to give an appropriate description of 
the pomeron. We might also remark that our calculation is done for zero 
transverse momentum Ax of the diffractively scattered proton, but a nonzero 
A|. would be small compared with the scale set by the p'^ of the final state 
jets. 

For heavy fiavour photoproduction the cancellation we have mentioned 
is not present and we find an appreciable cross section for charm quark jets 
in our model. If we assume that our result for light quarks is not radically 
changed by the corrections mentioned above, open charm should be domi- 
nating or at least enhanced in diffractive photoproduction. Provided that 
at least the part of the pomeron which couples to the jets is dominated by 
large gluon momenta, one might hope to study the hard pomeron in the case 
that experiment finds a significant rate of events with light quark jets, since 
they cannot be due to purely soft pomeron exchange. In any case it would 
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be interesting to study the flavour of the produced jets. 

From our calculation of gg-production we can estimate the total cross 
section for diffractive photon-proton scattering and extract from it the quark 
structure function of the pomeron. For this we need a particular ansatz for 
the gluon propagator D{P), but within a class of functions the results are 
approximately equal. For light quarks we flnd a structure function Gq/p{z) 
almost independent of Q^, similar in size to the one obtained in a more 
phenomenological approach to the pomeron (43) and with a shape somewhat 
harder than z[l — z). It is amusing to notice that in pp-coUisions, UA8 [20] 
have measured a pomeron structure function that is harder than zil — z), 
although with an uncertainty whether the dominating partons are quarks or 
gluons. Unlike in jet production, the dependence on the quark mass is rather 
strong, leading to a smaller structure function for strange quarks than for 
u and 0?, and to an uncertainty of a factor around 2 according to whether 
one takes current or constituent masses. For charm quarks the structure 
function is smaller by about an order of magnitude, and at large z it has a 
non-negligible contribution from longitudinal photons which is not scaling. 

We flnally note that some of our results have previously been obtained 
in the purely perturbative model of Nikolaev and Zakharov [9, 10, 11]. In 
particular, they have found that the integrated cross section for diffractive 
gg-production is suppressed for longitudinal photons and strongly flavour 
dependent if the photons are transverse, and that at large transverse quark 
momentum charm is enhanced with respect to light quarks. However, the 
relative contribution of charm is much smaller in their model than in ours, 
which was expected in [9]. Furthermore, the part of the quark structure 
function of the pomeron due to gg-production was found [10] to behave ap- 
proximately like zil — z). As for the production of light quark jets, it was 
reported [11] that the loop integration is dominated by squared gluon mo- 
menta of order p^ if the photon virtuality is very small compared with the 
invariant gg-mass, i. e. in particular if the photons are real. 
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A Results in a phenomenological model of 
the pomeron 

In this appendix we give the results of applying the phenomenological model 
of [13, 14] to the reaction investigated in the present paper. 

This approach to the pomeron builds on the observation that it behaves 
in many ways as an isoscalar photon, and its coupling to light quarks with 
small invariant mass is taken to be (3o^^. If one of the quark legs has large 
virtuality A;^, this coupling has to be multiplied by some form factor f{k^) to 
ensure that it vanishes as k^ — > — oo. The simple ansatz f{k^) = fi^ / {fi^ — k^) 
has been successful in describing the proton structure function at moderately 
low xb and diffractive p production in DIS. 

In this framework the two diagrams that describe the diffractive process 
p + 7 — > p + gg are those one obtains from fig. 1 when taking into account the 
two possible directions of charge fiow for the quark line. In the case of zero 
quark mass m^ equations (22) to (24) for the cross section then are replaced 
with 



dtd^dp^T 4^2 .'US L .v;. 3 + g2 /^31^ 



(/(t) + /(u))'x< 



1 — 2 —r for transverse photons 



s 



pi {s - Qr 



pI {s + 



12\2 



s Q'^s 



for longitudinal photons v / 



for scalar polarisation 



where by "scalar polarisation" we mean that the polarisation vector of the 
photon is set proportional to its momentum q. We remark in passing that 
this is very similar to what one would obtain in the LN model considering 
only diagrams (a) and (c) in fig. 3 (whose kinematics correspond to that of 
fig. 1).^"* The above result violates electromagnetic gauge invariance badly: 
the cross section for scalar photon polarisation is larger than that for longi- 
tudinal photons instead of being zero, and both are divergent when Q^ — > 0. 
For diffractive production of jets with large transverse momentum this model 
does not give physically reasonable results. Gauge invariance might be re- 
stored by introducing an appropriate contact interaction between the photon, 
the pomeron and the two quarks, as suggested by the field theoretic approach 
to diffractive scattering of ref . [21]. 



"'^■^The cross sections are identical apart from factors as(p|,)/as ' in the LN model and 
{f{t) -\- f{u)YlA in the form factor approach. 
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If one takes p'^ to be small compared with s, i. e. considers the case 
where the produced quarks are rather coUinear with the photon-pomeron 
axis, the longitudinal and gauge terms in (44) become small compared with 
the transverse one and we might believe the results to be more reliable. If 
we require that s -\- Q^ ^ fj,Q then 

and the integrated cross section for the process is dominated by the low p^ 
region (as it is in the LN model, cf. section 7). Integrating (44) over p'^ > we 
find that the results for longitudinal photons and pure gauge are suppressed 
compared with the transverse one by fJ,Q/{s + Q^) (apart from logarithms). 
From the leading term for transverse photons the quark structure function 
(43) of the pomeron has been extracted in [13]. 
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